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Abstract 

Background: Glioblastoma (GBM) is an aggressive disease associated with poor survival. It is essential to account 
for the complexity of GBM biology to improve diagnostic and therapeutic strategies. This complexity is best 
represented by the increasing amounts of profiling ("omics") data available due to advances in biotechnology. The 
challenge of integrating these vast genomic and proteomic data can be addressed by a comprehensive systems 
modeling approach. 

Methods: Here, we present an in silico model, where we simulate GBM tumor cells using genomic profiling data. 
We use this in silico tumor model to predict responses of cancer cells to targeted drugs. Initially, we probed the 
results from a recent hypothesis-independent, empirical study by Garnett and co-workers that analyzed the sensitivity 
of hundreds of profiled cancer cell lines to 130 different anticancer agents. We then used the tumor model to predict 
sensitivity of patient-derived GBM cell lines to different targeted therapeutic agents. 

Results: Among the drug-mutation associations reported in the Garnett study, our in silico model accurately predicted 
-85% of the associations. While testing the model in a prospective manner using simulations of patient-derived GBM 
cell lines, we compared our simulation predictions with experimental data using the same cells in vitro. This analysis 
yielded a -75% agreement of /n silico drug sensitivity with in vitro experimental findings. 

Conclusions: These results demonstrate a strong predictability of our simulation approach using the in silico tumor 
model presented here. Our ultimate goal is to use this model to stratify patients for clinical trials. By accurately 
predicting responses of cancer cells to targeted agents a priori, this in silico tumor model provides an innovative 
approach to personalizing therapy and promises to improve clinical management of cancer. 

Keywords: Glioblastoma, Cancer, In Silico modeling. Deterministic model. Virtual tumor technology. Tumor profiling. 
Personalized therapy. Targeted therapy 



Introduction 

Cancer remains a major unmet clinical need despite ad- 
vances in clinical medicine and cancer biology. Glioblastoma 
(GBM) is the most common type of primary adult brain 
cancer, characterized by infiltrative cellular proliferation, 
angiogenesis, resistance to apoptosis, and widespread gen- 
omic aberrations. GBM patients have poor prognosis, with a 
median survival of 15 months [1]. Molecular profiling and 



* Correspondence: skesari(5)ucsd.edu 
^Equal contributors 

^Translational Neuro-Oncology Laboratories, Moores Cancer Center, UC San 
Diego, La Jolla, CA 92093, USA 

^Department of Neurosciences, UC San Diego, La Jolla, CA 92093, USA 
Full list of author information is available at the end of the article 

(3 BioMed Central 



genome-wide analyses have revealed the remarkable gen- 
omic heterogeneity of GBM [2,3]. Based on tumor profiles, 
GBM has been classified into four distinct molecular sub- 
types [4]. However, even with existing molecular classifica- 
tions, the high intertumoral heterogeneity of GBM makes it 
difficult to predict drug responses a priori This is even 
more evident when trying to predict cellular responses to 
multiple signals following combination therapy. Our ration- 
ale is that a systems-driven computational approach will 
help decipher pathways and networks involved in treatment 
responsiveness and resistance. 

Though computational models are frequently used in 
biology to examine cellular phenomena, they are not 
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common in cancers, particularly brain cancers [5,6]. 
However, models have previously been used to estimate 
tumor infiltration following surgery [7] or changes in 
tumor density following chemotherapy in brain cancers 
[8]. More recently, brain tumor models have been used 
to determine the effects of conventional therapies in- 
cluding chemotherapy and radiation [5]. Brain tumors 
have also been studied using an agent-based modeling 
approach [9]. Multiscale models that integrate hierarch- 
ies in different scales are being developed for application 
in clinical settings [10]. Unfortunately, none of these models 
have been successfully translated into the clinic so far. It is 
clear that innovative models are required to translate data 
involving biological networks and genomics/proteomics into 
optimal therapeutic regimens. To this end, we present a de- 
terministic in silico tumor model that can accurately predict 
sensitivity of patient-derived tumor cells to various targeted 
agents. 

Methods 

Description of In Silico model (Version 7.3 Cellworks) 

We performed simulation experiments and analyses using 
the predictive tumor model - a comprehensive and dy- 
namic representation of signaling and metabolic pathways 
in the context of cancer physiology. This in silico model 
includes representation of important signaling pathways 
implicated in cancer such as growth factors such as EGFR, 
PDGFR, FGFR, c-MET, VEGFR and IGF-IR; cytokine and 
chemokines such as ILl, IL4, IL6, IL12, TNF; GPCR medi- 
ated signaling pathways; mTOR signaling; cell cycle 
regulations, tumor metabolism, oxidative and ER stress, 
representation of autophagy and proteosomal degradation, 
DNA damage repair, p53 signaling and apoptotic cascade. 
The current version of this model includes more than 
4,700 intracellular biological entities and -6,500 reactions 
representing their interactions, regulated by -25,000 
kinetic parameters. This comprises a comprehensive 
and extensive coverage of the kinome, transcriptome, 
proteome and metabolome. Currently, we have 142 kinases 
and 102 transcription factors modeled in the system. 

Model development 

We built the basic model by manually curating data from 
the literature and aggregating functional relationships be- 
tween proteins. The detailed procedure for model devel- 
opment is explained in Additional file 1 (Section 2) using 
the example of the epidermal growth factor receptor 
(EGFR) pathway block (Additional file 1: Figure SI and 
Figure S2). We have also presented examples of how 
the kinetic parameters are derived from experimental 
data, in Additional file 1: (Section 2). We have validated 
the simulation model prospectively and retrospectively, at 
phenotype and biomarker levels using extensive in vitro 
and in vivo studies [11-20]. 



Disease phenotype definitions 

Disease phenotype indices are defined in the tumor 
model as functions of biomarkers involved. Proliferation 
Index is an average function of the active CDK-Cyclin 
complexes that define cell cycle check-points and are 
critical for regulating overall tumor proliferation poten- 
tial. The biomarkers included in calculating this index 
are: CDK4-CCND1, CDK2-CCNE, CDK2-CCNA and 
CDKl-CCNBl. These biomarkers are weighted and their 
permutations provide an index definition that gives max- 
imum correlation with experimentally reported trend for 
cellular proliferation (based on literature). 

We also generate a Viability Index based on 2 sub- 
indices: Survival Index and Apoptosis Index. The bio- 
markers constituting the Survival Index include: AKTl, 
BCL2, MCLl, BIRC5, BIRC2 and XIAP. These biomarkers 
support tumor survival. The Apoptosis Index comprises: 
BAX, CASP3, NOXA and CASP8. The overall Viability 
Index of a cell is calculated as a ratio of Survival Index/ 
Apoptosis Index. The weightage of each biomarker is 
adjusted so as to achieve a maximum correlation with 
the experimental trends for the endpoints (based on 
literature). 

In order to correlate the results from experiments 
such as MTT Assay, which are a measure of metabolic- 
ally active cells, we have a 'Relative Growth' Index that is 
an average of the Survival and Proliferation Indices. 

The percent change seen in these indices following a 
therapeutic intervention helps assess the impact of that 
particular therapy on the tumor cell. A cell line in which 
the Proliferation/Viability Index decreases by <20% from 
the baseline is considered resistant to that particular 
therapy. 

Creation of cancer cell line and its variants 

To create a cancer-specific simulation model, we start 
with a representative non-transformed epithelial cell as 
control. This cell is triggered to transition into a neo- 
plastic state, with genetic perturbations like mutation 
and copy number variation (CNV) known for that spe- 
cific cancer model. We also created in silico variants for 
cancer cell lines, to test the effect of various mutations 
on drug responsiveness. We created these variants by 
adding or removing specific mutations from the cell line 
definition. For example, DU145 prostate cancer cells nor- 
mally have RBI deletion. To generate a variant of DU145 
with wild-type RBI (WT), we retained the rest of its muta- 
tion definition except for the RBI deletion, which was 
converted to WT RBI (Additional file 1). 

Simulation of drug effect 

To simulate the effect of a drug in the in silico tumor model, 
the targets and mechanisms of action of the drug are deter- 
mined from published literature. The drug concentration is 
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assumed to be post-ADME (Absorption, Distribution, Me- 
tabolism and Excretion). 

Creation of simulation avatars of patient-derived GBM cell 
lines 

To predict drug sensitivity in patient- derived GBM cell 
lines, we created simulation avatars {in silico profiles) for 
each cell line as illustrated in Figure IB. First, we simu- 
lated the network dynamics of GBM cells by using ex- 
perimentally determined expression data (Additional file 
1: Table SI; Additional file 1: Section 7). Next, we over- 
lay tumor-specific genetic perturbations on the control 
network, in order to dynamically generate the simulation 
avatar. For instance, the patient-derived cell line SK987 
is characterized by overexpression of AKTl, EGFR, IL6, and 
PI3K among other proteins and knockdown of CDKN2A, 
CDKN2B, RUNX3, etc. (Additional file 1: Table SI). After 
adding this information to the model, we further optimized 
the magnitude of the genetic perturbations, based on 
the responses of this simulation avatar to three mo- 
lecularly targeted agents: erlotinib, sorafenib and dasa- 
tinib. The response of the cells to these drugs (from 
in vitro experimental data) was used as an "alignment 
data set". In this manner, we used "alignment drugs" 
(erlotinib, sorafenib, and dasatinib) to optimize the 
magnitude of genetic perturbation in the trigger files 
and their impact on key pathways targeted by these 
drugs. For example, most GBM cell lines demonstrated 



dominance of EGFR signaling as they had gains in copy 
number of EGFR gene. Hence the effect of EGFR in- 
hibitor would be a good indicator for the relative dom- 
inance of this signaling pathway. This is illustrated in 
further details in Additional file 1 using an example of 
two cell line profiles that have EGFR over-expression 
but differential response to EGFR inhibitor. Similarly, so- 
rafenib helped determine and align with MEK/ERK activa- 
tion, while dasatinib with activation of SRC signaling. 

Simulation protocol 

The simulation protocol included 3 states: 

1. Control State - The in silico model was simulated 
for 50,000 seconds, during which the different 
biological entities (called species) attain a 
steady-state concentration. This concentration 
depends on the balance between the rate of 
reaction nodes producing the species and the 
reaction nodes utilizing/ degrading the species. 
This is an untriggered system and is representative of 
a non-transformed epithelial cell. 

2. Disease State - At 50,000 seconds simulation time, 
we introduced the mutation data (specific to 
patient-derived GBM cell lines to be created) and 
simulated for an additional 1,25,000 seconds. During 
this time, the system attained a new steady state that 
aligns to the network dynamics of the cell line. 
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Figure 1 Simulation worl^flow for in silico tumor model. A, This illustration depicts a representative simulation protocol used for retrospective 
analysis of gene mutation-drug sensitivity association reported in the Garnett study. The simulation protocol included 3 states: Control or 
Untriggered state is simulated for 50,000 seconds to allow the biological entities to attain a steady-state concentration. At 50,000 seconds, 
mutation data is introduced and simulated for an additional 1,25,000 seconds to attain Disease state. For Drug-treated state, we introduce a 
drug into the system by perturbing the target reaction nodes and simulate the model for 2,00,000 seconds. At the end of this state, we calculate 
percent change in the indices for cell survival. B, This schematic demonstrates the simulation workflow for creation, optimization and testing of 
patient-derived GBM cell line profiles in silico. The key steps involved in developing the simulation avatars of the patient-derived GBM cell lines include: 
Input profiling data reporting relative expression of the different proteins in the cell lines; Iterative testing and alignment of simulation avatars to match 
experimental data on drugs used for alignment of the network (erlotinib, sorafenib and dasatinib); Locking the simulation avatars; In silico predictions 
and /n vitro testing. 



Pingle et al. Journal of Translational Medicine 2014, 12:128 
http://www.translational-nnedicine.conn/content/1 2/1/1 28 



Page 4 of 1 3 



3. Disease State + Drug Treatment - Following the 
simulation run time of 1,25,000 seconds for Disease 
State, we introduced the drug into the system by 
perturbing the target reaction nodes as explained 
above. We then simulated the model further for 
2,00,000 seconds (drug treatment). A percent change 
in the indices for cell survival (described earlier) 
indicates the therapeutic potential of the drug. 
Iterative simulations with varying concentrations of 
the drug generate dose-response curves from which 
IC50 values can be determined. 

Figure lA is a schematic of the representative simula- 
tion protocol that we used for the retrospective analysis 
of gene mutations -drug effects reported in the study by 
Garnett and co-workers. Figure IB illustrates the work- 
flow for simulation studies on patient-derived GBM cell 
lines. For the patient-derived GBM cell line predictions, 
we prospectively compared in silico responses to experi- 
mentally obtained results {in vitro data from patient- 
derived GBM cell lines) and determined corroboration 
between in silico and in vitro data. As per the dose- 
response plots generated by in silico predictions, a cell 
line was considered sensitive to a drug if it demon- 
strated >20% decrease in relative growth. The 20% thresh- 
old was used for both in silico predictions and for in vitro 
experimental data. 

Patient-derived glioblastoma cell lines 

Fresh human glioblastoma samples were acquired from 
brain tumor patients undergoing clinically indicated sur- 
gery (University of California San Diego Human Subjects 
Protocol) and cultured as previously reported [21,22]. 
GBM4 and 8 cells were a kind gift from C. David James 
(University of California San Francisco). Briefly, the disso- 
ciated tissue was washed, filtered through a 30 (im mesh 
and plated onto ultra-low adherence flasks at a concentra- 
tion of 500,000 to 1,500,000 viable cells/ml. The stem cell 
isolation medium included human recombinant EGF 
(20 ng/ml), human bFGF (10 ng/ml) and heparin 
(2 (ig/ml). Sphere cultures were passaged by dissoci- 
ation using Acutase (Sigma), washed, resuspended in 
neural stem cell culture medium (#05750, Stemcell Tech- 
nologies), and plated on ultra low- adherence 96 well plates 
at 2000 cells per well for all subsequent drug testing. We 
characterized all patient-derived glioblastoma lines using 
histopathologic and integrated genomic analyses. The 
glioblastoma lines were profiled using the Aflymetrix 
Gene Chip Human Gene 1.0 ST Array. 

Drug screening 

Drug screens were performed on patient-derived GBM 
cell lines plated at 2000 cell per well in 96-well microtiter 
plates, incubated overnight. After 72 hours of incubation 



with drugs, cell viability was quantified by the Alamar Blue 
assay. Briefly, after incubation, Alamar Blue (#BUF012B, 
AbDSerotec) was added directly to the culture medium, 
and the fluorescence measured at 560/90 to determine the 
number of viable cells (Infinite M200, Tecan Group Ltd.). 

Results 

Our study involved a retrospective component where we 
predicted gene mutations - drug sensitivity associations 
defined in a recent hypothesis-independent study [23]. 
In addition, we predicted sensitivity of our profiled 
patient-derived GBM cell lines to targeted agents and 
compared these in silico predictions to in vitro experi- 
mental data. 

Retrospective validation of in Silico tumor model 

In the first part of the study, we evaluated the ability of 
the in silico tumor model to predict drug responses that 
were reported in the study by Garnett and colleagues 
[23]. A comparison of our predictions with the associa- 
tions reported in the Garnett study indicated the pre- 
dictive capability of our in silico tumor model. 

Our modeling library has definitions for 45 of the 639 
cell lines used in this study (Additional file 1: Table S2) 
and supports 70 of the 130 drugs studied (Additional file 1: 
Table S3). Further, we can represent 51 of the 84 genes 
screened for mutations (Additional file 1: Table S4). Of the 
448 significant gene mutation-drug response associations 
reported, our in silico model was able to accurately predict 
22 of the 25 testable associations from the Garnett study 
(>85% agreement; Additional file 1: Table S5). The gene 
mutation-drug response correlations from the Garnett 
study that are currently not supported by the system are 
listed in Additional file 1: Table S6. From the 25 gene mu- 
tation-drug response associations tested from the Garnett 
study (Additional file 1: Table S5), a few examples of the 
correlations are explained below. Figure lA depicts a rep- 
resentative schematic of this retrospective analysis using 
the simulation {in silico tumor model). 

BRAF Mutations and Drug Sensitivity 

The Garnett study showed that cells with BRAF mutation 
were sensitive to the MEKl/2 inhibitor AZD2644 [23]. To 
examine this association, we modeled cancer cell variants 
with wild-type BRAF in silico. Modeling data showed that 
cells with wild-type BRAF were resistant to AZD6244, 
when compared to the parent tumor cells with mutant 
BRAF. Thus, BRAF mutation conferred sensitivity to the 
MEKl/2 inhibitor in silico; this prediction validates the 
finding reported in the Garnett study (Figure lA). 40-60% 
melanoma patients carry BRAF mutations that activate 
MAPK signaling [24,25] and this association could have 
therapeutic implications for the treatment of patients with 
BRAF mutant melanoma. 
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Effect of different mutations on sensitivity to tyrosine 
Kinase inhibitors 

The Garnett study showed that cells with BRAF muta- 
tion were sensitive to the MEKl/2 inhibitor AZD2644 
[23]. To examine this association, we created cancer cell 
variants with wild-type BRAF in the in silico model 
Simulation data showed that cells with wild-type BRAF 
were resistant to AZD6244, when compared to cells with 
mutant BRAF. Thus, BRAF mutation conferred sensitivity 
to the MEKl/2 inhibitor; this validates the finding re- 
ported in the Garnett study (Figure 2A). 40-60% melan- 
oma patients carry BRAF mutations that activate MAPK 
signaling [24,25]. This association tested in Figure 2A may 
have therapeutic implications for the treatment of patients 
with BRAF mutant melanoma. 

ERBB2 (HER2) amplification is a biomarker for sensi- 
tivity to EGFR-family inhibitors [26]. In the in silico 



model, we tested for sensitivity to EGFR2 family inhibi- 
tors, lapatinib and BIBW2992. Specifically, we examined 
sensitivity of cancer cells in the presence of mutations 
and/or over-expression of BRAF, CDHl, ERBB2, CCNDl 
and MET. These predictions from simulations were com- 
pared with results obtained in the Garnett study and the 
predictive capability of our model was determined. 

In silico predictions indicate that BRAF mutation de- 
creases sensitivity of cells to lapatinib (Figure 2B), whereas 
CDHl mutant lines demonstrated higher sensitivity to 
lapatinib when compared to variants with wild-type 
CDHl (Figure 2C). Further, cMET over-expression showed 
increased sensitivity to lapatinib, as indicated by decrease 
in viability in cells with cMET over-expression (Figure 2D). 
Additionally, ERBB2 and CCNDl over-expression cor- 
related positively with lapatinib sensitivity (Additional 
file 1: Table S5). In all these simulation experiments 
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Figure 2 Retrospective analysis tests in silico predictions of gene mutations and sensitivity to EGFR family inhibitors. Associations 
reported in tlie Garnett study were tested in a blinded manner using our in silico model and predictions obtained were compared to results 
reported in the Garnett study. A, We created wild-type BRAF variants of four cancer cell lines - COLO205, HT29, MDAMB231 and U266 in silico 
and compared the effect of MEKl/2 inhibitor AZD2644 on these cell lines and on corresponding parent lines expressing mutant BRAF. Our data 
demonstrated that BRAF mutation increases sensitivity to AZD6244. B, We simulated three cell lines - HI 650, HI 975 and SW48 with wild-type or 
mutant BRAF and tested for sensitivity to the EGFR2 family kinase inhibitor, lapatinib. BRAF mutation decreases sensitivity of cells to lapatinib. 
C, Similarly, when four cell lines (AGS, H1437, MKNl and MKN45) were tested for sensitivity to lapatinib, we observed that CDHl mutation 
increases sensitivity to lapatinib. D, We generated cell lines with wild-type or MET over-expression and tested the effect of lapatinib (A549, 
AGS, H358 and HT29 cell lines). MET over-expression increases sensitivity to lapatinib. 
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testing sensitivity to lapatinib, our in silico predictions 
corroborated with associations reported in the Garnett 
study. 

CDKN2A mutation and drug sensitivity 

The Garnett study reported associations between tumor 
suppressor gene mutations and several anti-cancer drugs. 
We tested these associations in our in silico tumor model. 
In the in silico analysis, cells harboring wild-type CDKN2A 
were resistant to erlotinib whereas CDKN2A mutation was 
associated with erlotinib sensitivity (Figure 3A). Similarly, 
cell lines with mutant CDKN2A showed increased sensitiv- 
ity to dasatinib (Figure 3B), bortezomib (Figure 3C), and to 
the CDK4/6 inhibitor PD0332991 (Figure 3D). These pre- 
dictions/analyses from our simulation corroborated accur- 
ately with data from the Garnett study. 

Other gene mutation-drug response associations exam- 
ined in our simulation models are illustrated in Additional 
file 1: Table S5. In addition. Additional file 1: Table S6 lists 



correlations between gene mutations and drug responses 
reported in the Garnett study, which are currently not 
supported by our modeling technology. In spite of these 
limitations, we obtained -85% agreement of our simula- 
tion data with findings reported by Garnett [23]. 

Prospective evaluation of tumor model - patient-derived 
GBM cell lines 

Identifying drug sensitivities in tumors/cancers with different 
mutations is important for designing individualized therapies 
for cancer. To this end, we created in silico avatars of 8 
patient-derived GBM cell lines using genomic data (Methods 
and Additional file 1: Table SI) and predicted their sensitivity 
to various targeted therapeutic agents. We then tested these 
in silico predictions prospectively by comparing then with 
experimental data obtained by in vitro testing on the same 
patient-derived GBM cell lines (Figure IB). 

The patient-derived GBM cell lines were obtained from 
patient tumors resected surgically and cultured in vitro 
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Figure 3 Retrospective analysis evaluates CDKN2A mutation - drug response association by in silico modeling. Using simulation 
modeling, we tested the role of the tumor-suppressor protein CDKN2A on sensitivity to different inhibitors and compared these predictions to 
those reported in the Garnett study. A, Cells expressing mutant CDKN2A and their wild-type variants were simulated in the in silico tumor model 
for four lines - BxPC3, H1437, H1650 and SW48. CDKN2A mutation increased sensitivity of cells to erlotinib when compared to wild-type CDKN2A. 
B, Cells with mutant CDKN2A were more sensitive to dasatinib than cells with wild-type CDKN2A (A549, BxPC3, HCTl 16 and H460). C, COLO205, 
HT29, HI 437 and SW48 cell lines with mutant CDKN2A were sensitive to bortezomib more than cells expressing wild-type variants. D, CDKN2A 
mutant cells BxPC3, H1437, H1975 and HT29 also showed higher sensitivity to CDK4-Cyclin Dl inhibitor PD0332991 over the CDKN2A W variants. 
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(details in Methods). We have profiled these lines using 
Affymetrix Gene Chip Human Gene 1.0 ST Array. Using 
whole-exome sequencing, we recently tested the validity 
of these cells (maintained in cultures) for development 
and testing of personalized targeted therapies, based on 
their accurate representation of the original tumor profiles 
[27]. We have designated the different patient- derived 



GBM cell lines as: GBM4, GBM8, SK102, SK262, SK429, 
SK748, SK987 and SK1035. 

After generating in silico profiles of these cells, we opti- 
mized these simulation avatars in terms of strength of func- 
tional effect of the mutation on key pathways such as EGFR, 
RAS and Src/PBK. The rationale for this optimization is 
that expression data on these cells does not provide an 
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Figure 4 In silico modeling analysis and experimental in vitro data for drug responsiveness to 3 alignment drugs. A, Predictive dose 
response data for eriotinib witli percent cliange in viability. Cells showing decrease in viability of 20% or greater are considered sensitive to the 
drug. B, In vitro experimental results for effect of 1 eriotinib on viability in patient-derived GBM cell lines; viability was determined at 72 h 
using Alamar Blue assay. C, D, Predictive and experimental data for sorafenib. E, F, Predictive and experimental data for dasatinib. All drugs were 
tested in vitro at 1 pM. Dose-response curves for in silico data demonstrate the effects of increasing concentrations of the drugs - eriotinib, 
sorafenib and sunitinib on the viability of profiled patient-derived GBM cell lines in the simulation model. 
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Figure 5 (See legend on next page.) 
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(See figure on previous page.) 

Figure 5 In silico modeling and experimental in vitro data for drug responsiveness to tyrosine kinase inhibitors. Tiiis figure 
demonstrates in silico predictions of sensitivity and in vitro viability (respectively) in response to treatment with tyrosine kinase inhibitors: 
A, B, lapatinib, C, D, nilotinib, E, F, Imatinib and G, H, Sunitinib. Cells were exposed in vitro to 1 |jM tyrosine kinase inhibitors for 72 h and 
viability determined using Alamar Blue assay. The dose-response for in silico predictions is generated by iterative simulations with increasing 
concentrations of the drug in the model and the viability index is calculated. Cells showing decrease in viability of 20% or greater are considered 
sensitive to the drug. 



accurate measure of the dominance of different intracellular 
pathways. In order to interrogate this information on the 
pathways that play a dominant role in each tumor line (such 
as EGFR, RAS, PI3K, etc.), we used 3 anti-cancer agents (er- 
lotinib, sorafenib and dasatinib) targeting these pathways. 
This will achieve "alignment" and train the simulation ava- 
tars for further analyses (details in Additional file 1). The 
alignment for these 3 drugs could be best achieved in the 
following cell lines: GBM8, SK262, SK429, SK748, and 
SK1035. In cell lines GBM4 and SK987, there was a mis- 
match for sorafenib where the predictive trends were re- 
versed. GBM4 was sensitive to sorafenib experimentally but 
our in silico predictions showed it to be resistant; SK987 
was resistant experimentally but sensitive in predictive re- 
sults. Similarly, the experimental trend for SK102 resistance 
to dasatinib could not be met predictively. Correlation of 
predictive trends with alignment drugs is shown in Figure 4 
A-F. 

Predictions obtained by simulation modeling are pre- 
sented as dose-response plots for viability; decrease in 
viability of >20% was considered as sensitive. Experimen- 
tally, viability was determined by Alamar blue assay, in 
response to 1 [iM concentration of respective inhibitors 
at 72 h. These data represent viability as mean values 
from triplicate samples. 

We tested ten anti-cancer drugs in silico on the simu- 
lation avatars of the 8 patient-derived GEM cell lines in 
a blinded prospective study. These simulations generated 
predictions that we compared with in vitro experimental 
data (Additional file 1: Table S7A-D). Of the 80 in silico 
predictions, 61 (76.25%) predictions showed agreement 
with in vitro experimental results. Analysis of drug 
sensitivity correlation for all 8 GBM patient-derived 
cell lines, for all the 13 drugs is summarized in Additional 
file 1: Table S7. Figures 5A-H and 6A-H show a drug- 
wise comparison of in silico predictions (dose-response 
curves) and in vitro experimental results generated 
with testing 1 (iM concentration of each drug on these 
cell lines. 

Effect of tyrosine kinase inhibitors on patient-derived 
GBM cells 

For the EGFR family inhibitor lapatinib, simulation stud- 
ies predicted SK429, SK748 and SK1035 to be resistant, 
which were confirmed by in vitro data. Similarly, model- 
ing predicted GBM8, SK102, SK262 and SK987 to be 



sensitive and these predictions were in agreement with 
experimental data (Figure 5A and B). However, modeling 
predicted GBM4 to be resistant to lapatinib while in vitro 
data showed GBM4 to be highly sensitive to lapatinib 
(Figure 5B). For the tyrosine kinase inhibitor nilotinib, 
the model predicted GBM8 to be sensitive while all the 
other profiles to be resistant (Figure 5C). In vitro stud- 
ies demonstrated that GBM8 was indeed sensitive to 
nilotinib as predicted, but there was a mismatch with 
the experimental results for two lines - SK262 and 
SK1035. Experimentally, SK262 was found to be sensi- 
tive, whereas SK1035 was on the borderline of sensitiv- 
ity and resistance (Figure 5D). For imatinib, simulation 
predicted that all GBM lines except GBM8 were resist- 
ant (Figure 5E). The experimental results corroborated 
with this in silico prediction (Figure 5F). Sunitinib was 
the other multi-tyrosine kinase inhibitor tested. Our 
simulation predicted GBM8, SK102 and SK987 to be 
sensitive to sunitinib; however, only GBM8 was found 
to be sensitive in vitro, SK262 was predicted to be re- 
sistant to sunitinib but in vitro data found it to be moder- 
ately sensitive. On the other hand, GBM4, SK429, SK748 
and SK1035 were found to be resistant in both simulation 
and experimental data (Figure 5G-H). 

Effect of other drugs on patient-derived GBM cells 

Besides the tyrosine kinase inhibitors, correlation be- 
tween in silico predictions and experimental results for 
the 8 patient-derived GBM cell lines was also tested for 
drugs such as pitavastatin (HMG CoA reductase inhibi- 
tor), everolimus (mTOR inhibitor), celecoxib (COX2 in- 
hibitor) and bortezomib (proteasome inhibitor) (Figure 6 
A-H). For bortezomib, all profiles were predicted to be 
sensitive and these predictions matched with in vitro ex- 
perimental results (Figure 6A and B). For everolimus, 
in vitro results were in agreement with simulation pre- 
dictions for all lines except SK429 (Figure 6C and D). 
Our in silico model predicted GBM4, SK262, SK429, 
SK748 and SK1035 to be resistant to celecoxib; these pre- 
dictions matched with in vitro results. However, GBM8, 
SKI 02 and SK987 were predicted to show moderate sensi- 
tivity to celecoxib, but were found to be resistant in vitro 
(Figure 6E and F). For pitavastatin, the simulation pre- 
dicted 5 patient-derived GBM cell lines to be sensitive 
(GBM8, GBM4, SK102, SK262 and SK987), of which 
SK987 was found to be resistant in vitro. On the other 
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Figure 6 In silico modeling and experimental in vitro data for drug responsiveness to different drugs. Tiiis figure demonstrates in silico 
predictions of sensitivity and in vitro viability in response to treatment of patient-derived GBM cell lines with A, B, bortezomib, C, D, everolimus, 
E, F, celecoxib, and G, H, pitavastatin. All drugs were tested in vitro at 1 for 72 h and viability was assayed using Alamar Blue assay. Cells 
showing decrease in viability of 20% or greater are considered sensitive to the drug. 

V J 



hand, of the cell lines predicted to be resistant (SK429, 
SK748 and SK1035), SK1035 was sensitive in vitro and did 
not match with the prediction (Figure 6G and H), 

These data demonstrate a 76.25% agreement between 
in silico predictions of drug response and in vitro experi- 
mental data in patient-derived GBM cell lines. 

Discussion 

Developing an in silico model that takes into account 
the complex genotypes/phenotypes of cancer to accur- 
ately predict drug response will help personalize therapy 
with more efficiency. In this study, we developed and 
validated a virtual tumor model by retrospectively test- 
ing it against a dataset from a recent screening study 
[23]; we obtained a corroboration of -85% between our 
predictions and the results from this study. Following 
this retrospective validation, we generated in silico pre- 
dictions to prospectively test the sensitivity of patient- 
derived GBM cell lines to targeted agents. These analyses 
also demonstrated a high degree of agreement (>75%) be- 
tween in vitro experimental findings and in silico predic- 
tions. These studies validate our in silico tumor model 
and the simulation-based approach and provide critical 
proof-of-concept of a priori prediction of responses to tar- 
geted therapies. Thus, this model provides an effective 
platform for testing and developing personalized thera- 
peutic regimens for cancer patients. 

The genomic inputs that we used to create simulation 
avatars for patient-derived GBM cell lines were copy 
number variation data. A more comprehensive and ac- 
curate profile would require additional data (gene muta- 
tions, methylation status etc. along with copy number 
variation); this would help us develop a more representa- 
tive avatar and would likely improve the accuracy of our 
drug response predictions and provide higher correlation 
with experimental data. 

Genotypes of cancer cell lines have traditionally been 
used to correlate with drug sensitivity [28,29]. A similar 
recent study makes efficient use of gene expression pro- 
files to categorize colorectal cancers into different mo- 
lecular and clinically actionable subtypes [30]. Moreover, 
it is clear that using molecular tumor profiles to stratify 
patients for therapy affects response and progression- 
free survival [31]. However, increasing amounts of data 
from genomic, proteomic, transcriptomic and metabolo- 
mic profiling will likely require integration of these var- 
ied datasets and development of predictive systems 



modeling, which may hold the key to effective cancer 
therapy. 

Rapid screening of patient samples in real time with 
models such as the one we have developed can drive crit- 
ical therapeutic decision-making. Although our current 
model makes only cell-intrinsic predictions, we have been 
able to achieve a high rate of agreement between in silico 
predictions and in vitro findings. Future versions of this 
model are being refined to incorporate tumor microenvir- 
onment including aspects of angiogenesis, hypoxia, and 
tumor-associated inflammation. We believe that incorpor- 
ating these features into our model would more accurately 
represent the tumor in a patient. Importantly, this will fur- 
ther help improve our predictions for designing thera- 
peutic regimens for GBM patients. This model can also be 
adapted to identify potential mechanisms of resistance a 
priori and to design rational drug combinations that pre- 
vent emergence of resistance and development of escape 
pathways. 

Our in silico model aligns with NCI guidelines that 
emphasize evaluation of similar predictor models to de- 
termine their accuracy [12,15,32,33]. We intend to test 
this model in clinical trials and utilize it as a tool to ex- 
pedite clinical decision-making and determine drugs/ 
combinations most likely to benefit a patient. Addition- 
ally, models such as these will play important roles in 
testing new biological hypotheses. This is critical to the 
discovery of molecular drivers and critical networks in 
cancer pathophysiology and the development of better 
diagnostics and effective therapeutics. 
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